On the Origin of Zigzag Magnetic Order in Iridium Oxide Na2lrOs 
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We explore the phase diagram of spin-orbit Mott insulators on a honeycomb lattice, within the 
Kitaev-Heisenberg model extended to its full parameter space. Zigzag-type magnetic order is found 
to occupy a large part of the phase diagram of the model, and its physical origin is explained as due 
to interorbital ti g — e g hopping. Magnetic susceptibility and spin wave spectra are calculated and 
compared to the experimental data, obtaining thereby the spin coupling constants in Na2ir03 and 
Li 2 Ir0 3 . 
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In the quest for the materials with novel electronic 
phases, iridium oxide Na2lr03 came into focus re- 
cently (Tj-0] due to theoretical predictions [1, that this 
system may host Kitaev model physics and quantum spin 
Hall effect. 

Na 2 Ir03 is an insulator with sizable and temperature 
independent optical gap ~ 0.35 eV Q, and shows Curie- 
Weiss type susceptibility [l|,@ with moments correspond- 
ing to effective spin one-half of Ir 4+ ion with i|„ config- 
uration [lfj. These facts imply that Na 2 IrC>3 is a Mott 
insulator with well localized Ir-moments. 

Collective behavior of local moments in Mott insula- 
tors is governed by three distinct and often competing 
forces: (i) orbital-lattice (Jahn- Teller) coupling, (it) vir- 
tual hopping of electrons across the Mott gap resulting 
in exchange interactions, and (Hi) relativistic spin-orbit 



coupling (see Ref. for extensive discussions). The 



corresponding energy scales Ejt , J, and A vary broadly 
dependin g o n the type of magnetic ions and chemical 
bonding [12j . When A > (Ejt, J), as often realized 
for Co, Rh, Ir ions in octahedral environment, local 
moments acquire a large orbital component which may 
result in a strong departure from spin-only Heisenberg 
models @, [HI]. The direct observation of large spin- 
orbit splitting 3A/2 ~ 0.6 — 0.7 eV in insulating iridates 
Sr 2 Ir0 4 0, Sr 3 Ir 2 7 14], and Na 2 Ir0 3 made it 
certain that A > (Ejt, J)- Thus, low-energy physics of 
Na 2 Ir03 is governed by interactions among the spin-orbit 
entangled Kramers doublets of Ir-ions. 

It is also established now 0-01 that Ir-moments in 
Na 2 Ir03 undergo antiferromagnetic (AF) order at TV — 
15 K. The fact that Tjv is much smaller than param- 
agnetic Curie temperature (—125 K) Q and spin-wave 
energies 0] implies that the underlying interactions are 
strongly frustrated. This is natural in so-called Kitaev- 
Heisenberg (KH) model [l(| where long range order is 
suppressed by the proximity to the Kitaev spin-liquid 
(SL) state. However, the observed "zigzag" magnetic 
pattern [ferromagnetic (FM) zigzag chains, AF-coupled 
to each other] came as a surprising challenge to this sim- 
ple and attractive model. To resolve the "zigzag puzzle" , 



various proposals, ranging from routine extensions of the 
KH model 0, EE EH to a complete denial [H] of the 
presence of local Ir-moments, have been put forward. 

In this Letter, we show that the zigzag order is in fact 
a natural ground state (GS) of the KH model, in a previ- 
ously overlooked parameter range. Next, we identify the 
exchange process that supports a zigzag-phase regime. 
Further, we calculate spin-wave spectra and magnetic 
susceptibility of the model in zigzag phase, and find a 
nice agreement with experiments [J 0, @]. This lends 
strong support to the KH model as a dominant interac- 
tion in Na 2 Ir03 and related oxides. 

The model- Nearest-neighbor (NN) interaction be- 
tween isospin one-half Kramers doublets of Ir 4+ ions, 
coupled via 90°-exchange bonds, reads as follows (the 
exchange processes are described later): 



WW = 2K S]S] 



J S i ' $ j 



(i) 



Here, j(= x, y, z) labels 3 distinct types of NN bonds of 
a honeycomb lattice [16| of Ir ions in Na 2 Ir03 , and spin 
axes are oriented along the Ir-0 bonds of IrOg octahe- 
dron. The bond-dependent Ising coupling between the 
7 components of spins is nothing but Kitaev model [2(| , 
and the second term stands for the Heisenberg exchange. 

Let us introduce the energy scale A = V ' K 2 + J 2 and 
the angle tp via K = A sirup and J = A cos tp: the 
model |T|) takes then the following form: 



= A (2sin(^ S7S] + costp S v Sj). 



(2) 



We let the "phase" angle ip to vary from to 2-7T, un- 
covering thereby additional phases of the model that es- 
caped attention previously [16|, including its zigzag or- 
dered state which is of a particular interest here. 

It is highly instructive to introduce, following Refs. 11 
[l6j . 4 sublattices with the fictitious spins S, which are 
obtained from S by changing the sign of its two appro- 
priate components depending on the sublattice index. 
This transformation (generic for triangular, honeycomb, 
kagome lattices) results in the ^-Hamiltonian of the same 
form as ([T]), but with effective couplings K = K + J 
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FIG. 1: (color online), (a) Phase diagram of the Kitaev- 
Heisenberg model containing 2 spin-liquid and 4 spin-ordered 
phases. The transition points (open dots on tp-circle) are ob- 
tained by an exact diagonalization. The gray lines inside the 
circle connect the points related by the exact mapping (see 
text). Open/solid circles in the insets indicate up/down spins. 
The rectangular box in zigzag pattern (top-left) shows the 
magnetic unit cell, (b) Groundstate energy Egs of 24-site 
cluster and its second derivative — d 2 i?Gs/dv3 2 revealing the 
phase transitions. 



and J = — J, revealing a hidden SU(2) symmetry of the 
model at K = — J (where the Kitaev term K vanishes). 
For the angles, the mapping reads as tan <p = — tan ip — 1 . 

Phase diagram - In its full parameter space, the KH 
model accommodates 6 different phases, best visualized 
using the phase-angle ip as in Fig. JHa). In addition to 



the previously discussed [lQ, [2JJ, |22| Neel-AF, stripy- AF, 
and SL states near ip = 0, — j, and — ~, respectively, we 
observe 3 more states. First one is "AF" (K > 0) Ki- 
taev spin-liquid near tp= ~. Second, FM phase broadly 
extending over the third quadrant of the y-circle. The 
FM and stripy- AF states are connected [see Fig. QJa)] by 
the 4-sublatticc transformation, which implies their iden- 
tical dynamics. Finally, near tp = § 7t, the most wanted 
phase, zigzag-AF, appears occupying almost a quarter 
of the phase space. Thanks to the above mapping, it is 



understood that the zigzag and Neel states are isomor- 
phic, too. In particular, the <p = j7T zigzag is identical to 
Heisenberg-AF of the fictitious spins. 

To obtain the phase boundaries, we have diagonal- 
ized the model numerically, using a hexagonal 24-site 
cluster with periodic boundary conditions. The clus- 
ter is compatible with the above 4-sublattice transfor- 
mation and tp (p mapping. As seen in Fig. Hfb), 
the second derivative of the GS energy Eqs with re- 
spect to (p well detects the phase transitions. Three 
pairs of linked transition points are found: (87.7°, 92.2°) 
and (—76.1°, —108.2°) for the spin liquid/order transi- 
tions around ±5, and (161.7°, —33.8°) for the transitions 
between ordered phases. 

The transitions from zigzag-AF to FM, and from 
stripy-AF to Neel-AF are of first order by symmetry; 
see very sharp peaks in Fig. QJb). The spin liquid/order 
transitions near <p — — ^ lead to wider and much less 
pronounced peaks, suggesting a second (or weakly first) 
order transition [16j . On the contrary, liquid/order tran- 
sitions around tp = | show up as very narrow peaks; on 
the finite cluster studied, they correspond to real level 
crossings. Nature of these quantum phase transitions re- 
mains to be clarified. 

While at J = (i.e. tp = ±f ) the sign of K is irrel- 
evant [20] , the stability of the AF- and FM-type Kitaev 
spin-liquids against J-perturbation is very different: the 
SL phase near ^ (— f ) is less (more) robust. This phase 
behavior is related to a different nature of the competing 
ordered phases: for the ? SL, these are highly quan- 
tum zigzag and Neel states, while the SL near — § is 
sandwiched by more classical (FM and "fluctuation free" 
stripy [ll|) states which are energetically less favorable 
than quantum SL state. 

Exchange interactions in Na2lrO^,.- Having fixed the 
parameter space {K > 0, J < 0) for zigzag phase, we turn 
now to the physical processes behind the model (JXJ) . In- 
teractions between local moments in Mott insulators arise 
due to virtual hoppings of electrons. This may happen 
in many different ways, depending sensitively on chemi- 
cal bonding, intra-ionic electron structure, etc. The case 
of present interest (i.e., strong spin-orbit coupling, t\ g 
configuration, and 90°-bondi ng geom etry) has been ad- 
dressed in several papers d, [H], Hi H||. There are fol- 
lowing four physical processes that contribute to K and 
J couplings. 

Process 1: Direct hopping t' between NN t^g orbitals. 
Since no oxygen orbital is involved, 90°-bonding is irrel- 
evant; the resulting Hamiltonian is Hi — I\ Si-Sj with 
h ~ (|i') 2 /[/ 16]. Here, U is Coulomb repulsion be- 
tween tig electrons. Typically, one has t' '/t < 1, when 
compared to the indirect hopping t of t 2g orbitals via 
oxygen ions. 

Process 2: Interorbital NN t 2g — e g hopping t. This is 
the dominant pathway in 90°-bonding geometry since it 
involves strong tpda overlap between oxygen- 2p and e g or- 
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FIG. 2: (color online). Schematics of 4 different exchange pro- 
cesses (see text for details), arranged around the tp-phase dia- 
gram of Fig.[]Ja). Taken separately, the Hamiltonians Hi, H2, 
Hz, and H4 would favor "pure" Neel-AF, zigzag-AF, Kitaev- 
SL, and stripy- AF states, respectively, as indicated by arrows 
connecting Hi with the dots on ip-circle. The circle is divided 
into the phase-sectors by gray lines; SL phases are shaded. 



bitah typically, t/t ~ 2. The corresponding Hamiltonian 
is 



11]: 



H. 



(7) 



h (2SJS] 



Si ■ Sj). 



(3) 



This is nothing but the model (JTJ with K = — J = I2 > 0, 
i.e., at its SU(2) symmetric point ip = jtt inside the 
zigzag phase, see Fig. [2j For the Mott-insulating iridates 
(as opposed to charge-transfer cobaltates [Hj]), we es *i" 
mate I2 ~ |(i/(7) 2 Jh, where U is (optically active) ex- 
citation energy associated with t2 g — e g hopping, and Jh 
is Hund's interaction between ti g and e g orbitals. The 
physics behind this expression is clear: (t/U) 2 measures 
the amount of ti g spin which is transferred to NN e g or- 
bital; once arrived, it encounters the "host" ti g spin and 
has to obey the Hund's rule. 

For its remarkable properties, the Hamiltonian H2 (J3j> 
deserves a few more words. On a triangular lattice, it 
shows a nontrivial spin vortex ground state (see Fig. 5 of 
Ref. |11|): however, the elementary excitations are simple 
SU(2) magnons of a conventional Heisenberg-AF. When 
regarded as "J" -part of a doped t — J model, it leads to 
an exotic pairing [ll|, 0] ■ 

Process 3: Indirect hopping t between NN t2 g or- 
bitals via two intermediate oxygen ions. This gives 



rise to the Kitaev model H. 



(7) 



hS]S], with J 3 
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|(t 2 /J7)(Jff/J7) 8] where Jh is Hund's coupling between 
t2 g electrons. This process supports ip = — ? SL state, 
see Fig. EI 

Process 4 : Mechanisms involving pd charge-transfer 
excitations (energy A p( j). Two holes may meet at an 
oxygen (and experience U p repulsion), or cycle around 



a R2O2 plaquette (Fig. [2]). The corresponding Hamil- 
tonian H4 [8J, [ll|, [23j has the form as of H2 (0. The 
coupling constant I4 ~ |( A J^ v y 2 — -^-^) turns out to 
be negative because of the near cancellation of the two 
terms [1, [25| . It thus supports stripy- AF not observed in 
Na 2 Ir0 3 . 

Putting things together, we observe that it is the in- 
terorbital t 2g — e g hopping H 2 process that uniquely 
supports zigzag order. This implies also that multior- 
bital Hubbard-type models, when applied to iridates with 
90°-bonding geometry, must include e g states as well, 
even though the moments reside predominantly in the 
tig shell. 

Up to this point, we neglected trigonal field splitting 
A of the t2 g level due to the c-axis compression present 
in Na2lrC>3. This approximation is valid as long as A is 
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much smaller than spin-orbit coupling A ~ 0.4 eV 
I2H ] and seems to be justified, since the recent ab-initio 
calculations 19] suggest that A ~ 75 meV only 27|. 

We have also examined the longer-range couplings, 
using the hopping matrix of Ref. [19(, and found that 
second-NN interaction has the form of (J3j) (as noticed 
previously 28J), while third- NN interaction is of AF- 
Heisenberg type. The second (third)-NN interaction 
would oppose (support) zigzag order; however, we be- 
lieve that these couplings are not significant in Na2lrC>3 
because the corresponding long-range hoppings are found 

3 



to be small 

We do not attempt here to evaluate the parameters 
involved in H1-H4; ab-initio calculations as in Ref. [3(| 
might be more useful in this regard. Instead, having 
obtained a zigzag order in our model (JTJ) and identified 
the physical process driving this order, we turn now to 
the experimental data. The J and K values in Na2lrC>3 
and Li2lrC>3 will be extracted below from analysis of the 
neutron scattering and magnetic susceptibility data. 

Spin-waves in the zigzag phase - Consider a single do- 
main zigzag state, e.g., with FM chains running perpen- 
dicular to z-type bonds. Following Ref. 0], we introduce 
a rectangular a x b magnetic unit cell [\/3ao x 3ao in terms 
of hexagon-edge ao, see Fig.QTa)], and define the afr-plane 
wave vector q in units of (h, k) as q = (—h, nrfe)- Stan- 
dard spin-wave theory gives four dispersive branches: 

«? i2 (A, k) = [K 2 + (K + J) 2 ] c 2 - KJ(1 - s h s k ) 

±\(K + J)c h W{2K - J) 2 - (2Ks h - Js k ) 2 , (4) 



cosirh, Sh 

3„ 



and to^Qijk) = uii t 2(—h,k), with Ch 
sbnrh, and s/. = sin7rA:. If K = — J, i.e. at <p = |7r point 
of hidden SU(2) symmetry, two branches are degenerate 
(ui = L02) and become true Goldstone modes. Away 
from this special point, small magnon gap is expected 
to open by quantum effects not considered here. For q 
with h = k, the spin- wave dispersions Q simplify to 
wi(M) = yj2K{2K+J) \c h \ and u 2 (h,h) = V2\Jc h \, 
revealing two different energy scales in magnon spectra 
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(0,0) 



FIG. 3: Magnon spectra in the zigzag phase calculated using 
Eq. © with (J, K) = (-4.01, 10.45) meV. The inset shows the 
path along the symmetry directions in the reciprocal space; 
notations of Ref. 3 are used. 



set by K and J couplings. 

While the bandwidth of the lowest dispersive mode 
(set by J) is already known to be about 5-6 meV 
we are not aware of the high energy magnon data to 
estimate K in Na2lrC>3. We have therefore examined 
(see below) the magnetic susceptibility data [J Q , and 
obtained (J,K) = (-4.01,10.45) meV that well fit the 
susceptibility as well as neutron scattering data [3] • With 
this, we predict magnon spectra for Na2lrC>3 shown in 
Fig. [3] The lowest dispersive (J) mode is as observed 
indeed. However, mapping out entire magnon spectra is 
highly desirable to quantify the Kitaev term K directly. 

Magnetic susceptibility .- We have calculated the uni- 
form magnetic susceptibility x(T) of the model ((T|) on 8- 
site cluster (using full exact diagonalization) and 24-site 
cluster (using finite-temperature Lanczos method [3l|, 
l32l|). Both clusters are compatible with the zigzag or- 
der when periodic boundary conditions are applied. The 
parameters are varied such that J = Acostp is consis- 
tent with the neutron data Q while ip stays within the 
zigzag sector of Fig. Qla); this strongly narrows the pos- 
sible A"-window. For the data fits, we let ^-factor of Ir 4+ 
ion to deviate from 2 (due to the covalency effects 
and include T- independent Van Vleck term xo- The re- 
sult for J = -4.01 meV, K = 10.45 meV, g = 1.78, 
Xo = 0.16 x 10~ 3 cm 3 /mol fits the Na2lrC>3 data nicely 
(Fig.QJ; deviations occur at low temperatures only, when 
correlation length exceeds the size of the cluster used. 
The fit is quite robust: similar results can be found for 
small only variations, locating Na2lr03 near ip = 111±2° 
of the model phase diagram Fig [Ha). The spin couplings 
obtained are reasonable for the 90°-exchange bonds (as 
expected [Ulnl, they are much smaller than in 180°-bond 
perovskites [13. fl4|). The magnitude of Van Vleck term 
also agrees with our estimate %o — j^^b^a — 0.2 x 10~ 3 
cm 3 /mol for Ir 4+ ion, considering spin-orbit coupling 
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FIG. 4: Experimental magnetic susceptibilities for 

Na2lr03 (squares) and Li2lr03 [(| (circles) fitted by the- 
oretical results calculated with (J,K) = (—4.01, 10.45) meV 
and (—5.30, 7.85) meV, respectively. Exact \ °f the 8-site 
cluster is shown as solid lines. Lanczos results for the 24- 
site cluster are indicated by shading [3^]. Their comparison 
suggests that the calculated \ gives the thermodynamic limit 
down to T w 100 K where the finite-size effects become sig- 
nificant. 



A ~ 0.4 eV [H,|26|. 

For the sake of curiosity, we have also fitted x(T) 
data of Li 2 Ir03 @, a sister compound of Na 2 Ir03. Ac- 
ceptable results have been found for the angle window 
ip = 124 ± 6°; a representative plot for J = —5.30 meV, 
K = 7.85 meV, g = 1.94, xo = 0.14 x 10~ 3 cm 3 /mol is 
shown in Fig. [4j It is worth noticing that the value of 
J, which controls the bandwidth of the softest spin-wave 
mode (see Fig. appears to be similar in both com- 
pounds. This may explain why they undergo magnetic 
transition at similar TV ~ 15 K, despite very different 
high temperature susceptibilities. 

To conclude, we have clarified the origin of zigzag 
magnetic order in Na2lrC>3 in terms of nearest-neighbor 
Kitaev-Heisenberg model for localized Ir-moments. The 
model well agrees with the low-energy magnon and high 
temperature magnetic susceptibility data. A general im- 
plication of this work is that the interactions considered 
here should hold a key for understanding the magnetism 
of a broad class of spin-orbit Mott insulators with 90°- 
exchange bonding geometry, including triangular, honey- 
comb, hyperkagome lattice iridates. 
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